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Abstract — Surface reconstruction is a very challenging step in 
Reverse Engineering. It generates a surface from point cloud 
acquired from a part surface. NURBS surfaces are commonly 
used for freeform surface reconstruction. There are several 
algorithms of NURBS surface reconstruction including: Least 
Squares Methods, simulated annealing, and particle swarm 
optimization. The major problem of these methods is that they 
require parameters and knots optimization which increases the 
complexity of the problem. The purpose of this paper is to 
introduce basis function neural networks which eliminates the 
need for parameters or knots optimization besides providing 
acceptable error. They are called Rational B-spline Neural 
Networks (RBNN). Also, they have the advantage of providing the 
control points and weights of approximated NURBS surface over 
regular function approximation networks. Training of the 
network is done using the back-propagation algorithm. Results 
showed that the RBNN have better approximation to NURBS with 
acceptable error provided that the number of control points and 
training rate are selected properly. 

Keywords — Freeform Surface, Reverse Engineering, Surface 
Reconstruction, NURBS, Neural Networks 



by interpolation or approximation as shown in fig.l. In 
interpolation, the resulting surface passes through the point 



I. 



Introduction 



Reverse Engineering (RE) is the process of transforming 
real parts into digital geometric models. These models can be 
then modified, analysed, and used for manufacturing these 
parts using CAD/CAM applications. RE is realized in two steps: 
surface digitizing and surface reconstruction. Surface digitizing 
measures the part surface using contact or noncontact 
techniques. It provides point cloud describing the surface 
geometry [1,2,3,4]. 

Surface reconstruction can be defined as the process of 
recovering the 3D shape of a part. It can be done by: polygonal 
mesh, or CSG model, or surface fitting. The models produced 
by polygon mesh are not suitable for CAD/CAM applications 
as they are not in an engineering drawing format. They are 
used for rapid prototyping, and animations [1, 5]. Also, the 
CSG model is suitable only for prismatic shapes and cannot 
deal with freeform surfaces [15, 21]. 

Surface fitting is the process of constructing a surface that 
has the best fit to a series of points. Surface fitting can be done 







Interpolation 




Approximation 
Fig. 1 Difference between Interpolation and Approximation 

cloud. This is useful when measurement is done at high 
precision. In approximation, the produced surface passes as 
closely as possible to the point cloud. This means 
approximation is more suitable for point clouds to smooth 
residual noise presented in the points and to reduce the 
computational effort required by interpolation. 

Parametric surfaces are commonly used for surface fitting. 
They include Bezier B-Spline and NURBS. NURBS surfaces 
are preferred in engineering design and manufacturing 
applications because they permit simple object-shape 
modification by changing only a small number of parameters, 
such as control points, knot vectors, and weights [6,7,13]. 



JCET Vol.l No.l July 2011 PP.34-38 www.ijcet.org ©World Academic Publishing 



Journal of Control Engineering and Technology (JCET) 



In this context, the paper presents a new method for 
NURBS surface approximation without the need of parameters 
or knots optimization. This method is based on the use of basis 
function neural networks. They are called Rational B-spline 
Neural Network (RBNN). They are used due to their higher 
approximation ability to NURBS surface and they provide the 
control points and weights which cannot be done by regular 
function approximation networks. Also, an algorithm is 
presented for training RBNN based on the back-propagation 
training algorithm. 

II. Previous Work 

Supposing that there is a point cloud containing (r x s) {Qik} 
points and it is required to approximate these points to a 
NURBS surface SCu/ jfo Vj_/J. Assuming that S(ui k ,v t>k ) has (p ,q) 
degrees with (m x n) control points in u and v directions, it can 
be defined as: 
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where {Py} are control points forming the control grid, 
{wy} are the weights, {N ip (u)} and {JV J? (VJ}are the non-rational 
B-spline basis functions defined over open non uniform 
knot vectors: 
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where r = n+p+1 and s = m+q+1. 

NURBS approximation means determining the control 
points and weights of S(ui >k ,VijJ which minimize the distance 
between the original points and their corresponding points on 
S(uijoVijt). This can be formulated as a least square expression 
as following: 

E=E[ =1 ZU(Qi,k-S0w ()k )) 2 

Before applying approximation methods, a 
parameterization technique should be used to assign parameter 
values (U( j foV;,( r ) to each point. Chord length method is the most 
common method. It can be formulated as following: 
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After parameterization, knot vectors are to be determined. 
Knot vectors are used with parameters to calculate the B- 

spline basis function. To reflect the distribution of parameters, 
the averaging method is employed to determine the knot 
vectors. The internal knots can be calculated by: 
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Once parameters and knot vectors are calculated, an 
approximation method is applied to compute the control points 
and weights of S(ui k ,VijJ . 

Approximation methods for NURBS surface are divided 
into two categories: direct methods, and surface skinning 
methods. Direct methods consist of simultaneously estimating 
the unknown control points of the surface and weights. Surface 
skinning, also known as lofting, is a process of passing a 
smooth surface through a set of so called cross-sectional curves. 
These curves approximate the measured points. They must 
have the same degree and knot vectors. This leads to high 
number of control points making the process more complex 
and time consuming than direct methods [9,10]. 

Direct methods include: Least Square Methods (LSM), 
optimization techniques, and neural networks. The LSM is still 
the most commonly used method. Weiyin Ma and J P Kruth 
[11] used the LSM to approximate measured points to a B- 
Spline surface. Ma and Kruth [12] used two LSM for NURBS 
curves and surfaces fitting which automatically identify the 
control points and their respective weights. The weights of 
control points are first identified through singular value 
decomposition, the control points are then determined by least 
squares minimization. It was found that the surface error 
produced by LSM was higher. Djordje Brujic et. al. [13] used 
LSM modified by sparsity structures of the relevant matrices to 
the problem and regularization terms. This improves the 
computational complexity and saves time producing lower 
errors. 

As parameterization and knot determination are done by 
approximated methods, these are reflected in the surface error. 
To overcome this problem, work has been done in two 
directions. The first is developing more accurate methods. 
Aziguli Wulamu et. al. [14] developed a parameterization 
called "Adaptive Equidistant Parameterization". The surface 
errors obtained from this algorithm was lower than those of 
chord length and equidistant methods. Kunal Soni, Daniel 
Chen, and Terence Lerch [15] used feature based 

parameterization. 

The second direction was to calculate initial parameters and 
knots by conventional methods, then optimization techniques 
were then used to refine them at each iteration. Jiing-Yih Lai, 
Chiou-Yuan Lu [5] used the cumulative chord-length method 
to assign the initial parameters to the points. Then the Powel's 
method is applied to optimize the parameters at each iteration 
until the desired error is reached. Mohammed Riyazuddin [16] 
used LSM to find the control points of the NURBS surface. 
Then simulated annealing was applied to optimize the values of 
weights then the knot vectors. 

Meerja Huma [17] used the Simulated Evolution algorithm 
to optimize the values of knots and weights. Akemi Galvez et. 
al. [18] used genetic algorithm to find the parameters of the 
points. Akemi Galveza et. al. [19] used two step genetic 
algorithm. The first is for parameterization and the second for 
knots determination. Nallig Leal et. al. [20] used evolutionary 
strategy to obtain the weights of the NURBS surface so that the 
distance between the point cloud and the NURBS is minimal. 
Akemi Galvez and Andres Iglesias [21] used particle swarm 
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Fig. 2 Architecture of surface approximation RBNN 



optimization with LU decomposition to compute parameters, 
knot vectors, control points, and weights. 

Neural networks are used for surface approximation due to 
their ability to function approximation. Regular function 
approximation neural networks include: Multi-Layer Feed- 
Forward (MLF), and Radial Basis Function (RBF). P Gu and X 
Yan [22] used MLF neural network to approximate 
parameterized points generated from NURBS surface. Ju H. et. 
al [4] used RBF neural network to approximate parameterized 
points acquired by laser scanning freeform surface. After 
training, these networks can be used only to calculate points on 
the approximated surface. Although MLF and RBF neural 
networks have higher approximation ability, they do not 
provide a solution capable of integration with CAD/CAM 
applications. This is because they do not provide the control 
points and weights of the surface which define the parametric 
surface. Also, using the generated points for data transfer 
between packages is more difficult and slower than using the 
parameters of the surface. 

To overcome the stated problems, parametric basis function 
neural networks were developed. They are similar in 
architecture to MLF neural networks. They have fast 
approximation ability to approximate points to parametric 
curves or surfaces due to using parametric basis functions as 
activation functions. Also, the NURBS control points and 
weights can be provided after training to reconstruct the 
surface in CAD/CAM applications, as they are represented in 
the network by the connecting weights. George K. Knopfa and 
Jonathan Kofman [7] first used this type of neural networks to 
approximate points to Bezier surface. Xiaogang G. [8] used B- 
Spline basis function neural network to approximate points to 
B-Spline curves and surfaces. 

III. RBNN 

RBNN are modified B-Spline basis function neural 
networks which can be used to approximate points to NURBS 
curves or surfaces. It is noted that the special architectures 
enable the network inputs and outputs to have the same 
relationship between the NURBS parameters and coordinates. 



There are two architectures of RBNN: the first is for NURBS 
curve approximation and the second is for NURBS 
surface approximation. 

A. Surface RBNN Architecture 

The architecture of surface approximation RBNN is shown 
in Fig. 2. This architecture has two inputs, parameters u lk and 
Vik of a point, and three outputs which are approximated Xcyo 
Ycik, and Zc lk coordinates of the point. The network has five 
layers. The number of neurons in the first layer equals to the 
sum of the number of neurons (m + n) in both directions u and 
v. The second and third layer have (m x n) neurons, the forth 
layer has four neurons, and the fifth layer has 3 neurons. The 
net input functions are linear sum in the first, third, and forth 
layers while the second and fifth have delta function. The 
activation functions are linear for all layers except the first has 
B-Spline basis function. The connecting weights between the 
second and third layers represent the weights of the NURBS, 
Wy, and they are updated at every training cycle (epoch). The 
connecting weights between neurons in the third layer and the 
first, second, and forth neurons in the forth layer represent the 
X, Y, and Z coordinates of the control points (Pxj j - Pyy - Pz; j) 
respectively. These connections are also updated at every 
epoch. The rest of connections equal unity and are fixed during 
training. By starting the networks, it performs linear sum of the 
terms of the networks from the inputs to the outputs. The 
outputs are calculated by: 
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B. Training 

During training, the connecting weights of the network are 
iteratively adjusted to minimize the network performance 
function. The performance function for the network is the mean 
square error between the network outputs and the target outputs. 
The outputs are the approximated Xc lk , Yc^k, and Zci_k 
coordinates of the points, while the target outputs are the 
measured Xt 1?k , Ytyo and Zti ( k coordinates of the points. For (r 
x s) points, the performance function can be calculated by: 

r s 

E = 2 2 ((Xt kk - Xc lik f + (Yt lik - Yc lik f 

1=1 k=l 

+ (Zti, k -Zc lk ) ) 
By applying the back-propagation training algorithm, the 
training rules for epoch k+1 can be expressed as: 
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where T| represents the training rate. The training rate is held 
constant throughout the training. The performance of the 
algorithm is very sensitive to the proper setting of the training 
rate. If the training rate is set too high, the algorithm may 
oscillate and become unstable. If the training rate is too low, 
the algorithm will take too long to converge. To obtain a better 
training process, an algorithm was developed. Fig. 3 shows the 
flow chart of the algorithm followed for training the network to 
reach the desired error. This algorithm starts by using small 
training rate and number of control points. Then they are 
gradually increased at each epoch until reaching the desired 
error. 

IV. Case Study 

To validate the developed algorithm, it was applied on 600 
points, shown in Fig. 5-a, obtained from B-spline surface. 
These points were to be fitted to a bicubic NURBS surface. At 
first, the network was constructed and the initial number of 
control points was assumed to be (4x4). The initial value of 
training rate was set to (0.01). By training the network, it was 
found that the mean square error was high. By increasing the 
number of control points and the training rate, the mean square 
error begin to drop. When the number of control points was 
(11x11) and the training rate was (0.1169), a mean square error 
of (15.6 xl0" 3 mm) was reached. This error was considered as a 
suitable error according to the number of points and the large 
area of the surface form which points were obtained. The 
approximated NURBS surface is shown in Fig. 4-b. Plot of the 
training process is shown in Fig. 5. Training was done in 100 
epochs in 18 seconds. The code was written in MATLAB. 

V. Conclusion 



Freeform surface reconstruction is a major challenge in 
Reverse Engineering. The goal of surface reconstruction is to 
approximate point cloud acquired from the surface to a suitable 
surface model with acceptable error. NURBS surfaces are 
widely used in the approximation of point clouds for freeform 
surfaces. NURBS approximation methods are used to compute 
the control points and weights of the approximated surface. 
The problem of most of these methods is that they require 
parameters and knots optimization to provide more 
enhancements to reduce the approximation error. This is 
because the methods of calculating knot vectors and parameters 
are approximated methods. Optimization adds more 
complexity to the process and makes it more difficult to solve. 
In this paper, RBNN are used for NURBS surface 
approximation. Compared to other NURBS approximation 
methods, they do not require parameters or knots optimization. 
This is due to the high ability of function approximation of 
neural networks. Also, after training RBNN, they provide the 
control points and weights of the NURBS surface which is not 
available in regular function approximation networks. The 
performance of the network depends on two factors: the 
training rate and the number of control points. It is found that 
the optimum value of training rate depends on the volume of 
point clouds, and the number of control points. Large volume 
point clouds require higher training rate while smaller require 
smaller training rate. Also, increasing the number of control 
points should be accompanied with an increase in the training 
rate. Therefore, the values of the number of control points and 
training rate should be balanced together to achieve optimum 
performance for the network. 



Start 
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of control points 
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Fig. 3 Flow chart of the algorithm followed for training RBNN to reach 
the desired error 
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(b) 

Fig. 4 (a) 600 regular data point (b) Approximated NURBS surface. 
Training of the RBNN 
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Fig. 5 Training of the network 
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